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Abstract 



Assuming the existence of a heavy Majorana neutral particle arising from a com- 
posite model scenario we discuss the constraints imposed by present experimental 
limits of half-life neutrinoless double beta decay {OujSjS) measurements on the cou- 
pling of the heavy composite neutrinos to the gauge bosons. For neutrino masses 
Mjv = 1 TeV we obtain a rather weak lower bound on the compositeness scale: 
A > 0.23 TeV. 



Heavy neutral Majorana particles with masses in the TeV region are predicted 
in various theoretical models, such as superstring-inspired Eg grand unification [|I| 
or left-right symmetric models . In addition the possibility of a fourth generation 
with a heavy neutral lepton, that could be of Majorana type, is not yet ruled out 

ii- 

In this paper we discuss the possibility that a heavy Majorana neutrino might 
arise from a composite model of the ordinary fermions |^ . Composite models, which 
describe quarks and leptons as bound states of still more fundamental particles, 
generally called preons, have been developed as alternatives to overcome some of 
the theoretical problems of the standard model ^ . 

Although no completely consistent dynamical composite theory has been found to 
date, various models have been proposed, and one common, (inevitable), prediction 
of these models is the existence of excited states of the known quarks and leptons, 
much in the same way as the hydrogen atom has a series of higher energy levels above 
the ground state. The masses of the excited particles should not be much lower than 
the compositeness scale A, which is expected to be at least of the order of a TeV 
according to experimental constraints. For example the search for four-fermion 
contact interactions gives A{eell) > 0.9 — 4.7 TeV depending on the chirality of the 
coupling and on the lepton fiavour |^, ^. We expect therefore the heavy fermion 
masses to be of the order of a few hundred GeV. The CDF experiment has excluded 
excited quarks in the mass range 90 — 570 GeV from 7 + jet and W + jet final states 
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Phenomenological implications of heavy fermions have been discussed in the liter- 



ature ||T0|, [Tj using weak isospin (Iw) and hypercharge (Y) conservation. Assuming 
that such states are grouped in SU{2) x f/(l) multiplets, since light fermions have 
Iw = 0, 1/2 and electroweak gauge bosons have Iw = 0, 1, to lowest order in pertur- 
bation theory, only multiplets with Iw < 3/2 can be excited. Also, since none of the 
gauge fields carry hypercharge, a given excited multiplet can couple only to a light 



inultiplet with the same Y. In addition, current conservation forces the couphng of 
the heavy fermions to gauge bosons to be of the magnetic moment type. 
We will only consider here the excited multiplet with Iw = 1/2 Y = —1 

which can couple to the light left multiplet 

through the gauge fields W'^a.nd B'^, with the additional assumption that N is a 
neutral Majorana fermion. 

In terms of the physical gauge fields W^^ = {l/\/2)[W^ =F i^^) the relevant 
effective interaction can be expressed as 



^'" = (4r){(^'''"^-)«^^''. 



+ (^Ea^'' — z/) d^W~ + h.c] + neutral currents. (3) 

where / is a dimensionless coupling constant, A is the compositeness scale, and r 
are the Pauli SU{2) matrices, and the rest of the notation is as usual in the standard 
model. An extension to quarks and other multiplets, with a detailed discussion of 
the spectroscopy of the excited particles can be found in Ref. [1TT| . 

Regarding the experimental mass limits on the heavy Majorana neutrinos from 
pair production, Z -^ NN, we have Mj^ > 34.6 GeV at 95% c.l., which has been 
deduced from the Z line shape measurements |T^ , and which is independent of the 
decay modes. More stringent limits ~ 90 GeV come from single excited neutrino 
production, Z -^ Nu, through the transition magnetic coupling, but these do depend 
on assumptions regarding the branching ratio of the decay channel chosen [P, |9|, p!^ . 

In practical calculations of production cross sections and decay rates of excited 
states, it has been customary [jl5|, |l^, [l^ to assume that the dimensionless coupling 
/ in Eq. (3) is of order unity. However if we assume that the excited neutrino 



is of Majorana type, we have to verify that this choice is compatible with present 
experimental limits on neutrinoless double beta decay (Oz//?/3): 

A{Z) -^ A{Z + 2) + e- + e- (4) 

a nuclear decay, see Fig. 1, that has attracted much attention both from particle 
and nuclear physicists because of its potential to expose lepton number violation. 
More generally, it is expected to give interesting insights about certain gauge the- 
ory parameters such as leptonic charged mixing matrix, neutrino masses etc. The 
process in Eq.(4), which can only proceed via the exchange of a massive Majorana 



neutrino, has been experimentally searched for in a number of nuclear systems |]T8 



and has also been extensively studied from the theoretical side [19, ^, 21]. 
We will consider here the decay 

76Ge ^ 76gg ^ 2e- (5) 

for which we have from the Heidelberg-Moscow /3/?-experiment the recent limit p2 



(Ti/2 is the half life = log2 x lifetime) 

Ti/2 ( ^^Ge ^ ^^Se + 2e^) > 1.95 x 10^^ yr 90 % c.l. (6) 

In the following we estimate the constraint imposed by the above measurement 
on the coupling (//A) of the heavy composite neutrino, as given by Eq. (3). The 
fact that neutrinoless double beta decay measurements might constraint composite 
models, was also discussed in ref . § but within the framework of a particular model 
and referring to a heavy Majorana neutrino with the usual 7^ coupling. 

The transition amplitude of Oz//3/3 decay is calculated according to the interaction 
lagrangian: 

9 \fi 



■'int 



+ cos ecJl:{x)W^^-\x) + h.c.} (7) 



where 9c is the Cabibbo angle (cos 9c = 0.974 ) and Jjl is the hadronic weak charged 
current 

k 

j,{k) = 'N{vk)i,Uv-fAl^)r-{k)M{rk) (8) 

and where r^ is the coordinate of the k-th nucleon, M = it'') and t_(A;) = 
(l/2)(ri(A;) — iT2{k)) is the step down operator for the isotopic spin, (r(A;) is the 
matrice describing the isotopic spin of the k-th nucleon) . We emphasize that in Eq. 
(7) we have a a^u type of coupling as opposed to the 7^ coupling so far encountered 
in all Ou/SP decay calculations. 

For simplicity, we carry out our analysis assuming that there are no additional 
contributions to Oi//?/? decay from light Majorana neutrinos, right handed currents 
or other heavy Majorana neutrinos originating from another source. 

The transition amplitude is then 
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-j={l - Pl2MPl)(r^,x{l + 75) p_^^ (l + I5)(rupv{p2) x 

[F{Z + 2, e,)F{Z + 2, e2)f'e'^-^e'P-yfA{{k - Pif)fA{{k + p^f) x 

^ "^'^ ^ ^ "^'^ [{k - p.? - MIA \{k + P2r- M^] ^'^ 

where (1 — P\2)l \/2 is the antisymmetrization operator due to the production of two 



identical fermions, the functions F{^Z^ e) are the well known Fermi functions p3 



that describe the distorsion of the electron's plane wave due to the nuclear Coulomb 
field (ej are the electron's kinetic energies in units oirrieC^), 

FiZ,e) = xiZ,e)j-^^^^ (10) 

27iaZ 
e 

and the nucleon form factor, 



x{Z,e) ~ X ' '{Z) = -^ — — — (Rosen-Primakoff aproximation) 



(1 + |qr/mji)^ 
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with rriA = 0.85 GeV, is introduced to take into account the finite size of the nucleon, 
which is known to give important effects for the heavy neutrino case. 



As is standard in such calculations, we make the following approximations: [19, 



i) the hadronic matrix element is evaluated within the closure approximation 

< fmx) r,iy)\t > ^ e^(^/-<^"»-"e^«^">-^»)^« < /|4^(x) r,{y)\t > (12) 

where < En > is an average excitation energy of the intermediate states. This allows 
one to perform the integrations over fco, Xq, y^ in Eq. (9); 

ii) neglect the external momenta pi, p2 in the propagators and use the long wave- 
length approximation e"'^^'^ = e"*^^'^ ~ 1; 

Hi) the average virtual neutrino momentum < |k| > ^ I/-R0 = 40 MeV is much 
larger than the typical low-lying excitation energies, so that, ko = Ef + Ei— < En > 
can be neglected relative to k; 

iv) the effect of W and N propagators can be neglected since M^ ^ 80 GeV is much 
greater than |k| in the region where the integrand is large, and we are interested in 
heavy neutrino masses M^ 3> Mw- 



Using the same notation as in Ref. |20l] we arrive at 



% = {GFCOs9c)'^h7rS{Eo~E,-E,)x 

-^(1 - Pi2)M(pi)a,„(T,, (1 + ^5HP2) [F{Z + 2, ti)F{Z + 2, £2)]'^' x 
^^E^^.</l/(^)/(OK> (13) 

kl 

where lij is an integral over the virtual neutrino momentum, {rki = r^ — ri , ru = 
\vk-ri\,Xki = mArki) 



'' Ml J (27r)3 {l + \k\'/m\Y 



with: 



S^JFA{xkl) + '-^^^^Fb{xm) (14) 



Fa{x) 
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[x + x) 



—X ^3 



M^) = ^(^ "^ 



(15) 



Since lij is a symmetric tensor, we can make the replacement (y^ia^j -^ 
{l/2){a^i,auj} = i]f,i,r]ij - i]i^i]ii_, + i^^^^^iuj- Then, using the nonrelativistic hmit 
of the nuclear current 



3,{k) 



(16) 



fvr-{k) if/i = 

-fAT~{k){(yk)i iin = i 
((Tfc is the spin matrice of the k-th nucleon) we arrive, with straightforward algebra, 
at 



S 



fi 



Mf, 



Mf,2nS{Eo-Ei-E2) 



(17) 



where we have defined 
1 



/ 
< m > 

Vn 



T2 



1 - Pi2)u{pi){l + -i^)v{p2) [F{Z + 2, ei)F{Z + 2, es)] 



1/2 



ITT'cVN < / I ^ U > 



m, 



mpime i^^i 



Er4A:)r_(0- 
r-ki 



'I' 

J' A 



7j - o'a: ■ <?« {Fb{xm) - 3Fa{xm)) 



-Ok ■ oi Fa\xu) H 2 Fb\xu] 

'^kl 

and Rq = tqA^^^ is the nuclear radius (ro = 1.1 fm). 

The new result here is the nuclear operator Q which is substantially different 

from those so far encountered in 0z//9/9 decays, due to the ct^j, coupling of the heavy 

neutrino that we are considering. The decay width is obtained upon integration 

over the density of final states of the two-electron system 

final spins 



dV 



{2nf2Ei{27if2E2 



(19) 
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and the total decay rate F can be cast in the form 

r = {G,cosecr^-^^^i^foAeo,z)\nM' (20) 

fou = ^Oufou ' (21) 

foY- = \x''-''iZ + 2)\'^{et + 10el + AOel + 60eo + 30) (22) 

where, Qfi = < f\^\i >, eo is the kinetic energy of the two electrons in units of rrieC^, 
and ^Oi^ is a numerical factor that corrects for the Rosen-Primakoff approximation 
pO| used in deriving the analytical expression of f^^'- For the decay considered in 
Eq.(5), we have [^ C,o„ = 1.7 and eo = 4. The half-life is finally written as 



AV = (l0S2) ,^ ''f ,J^'f = 1.24 X 10'«yr 



Combining Eq. (23) with the experimental limit given for the decay considered 
in Eq. (5), we obtain a constraint on the dimensionless coupling / 



\f\< 



'M^A^y/' 



ni„ni 



p"''A, 



1/4 ^ 



fif,P/2 



'/i 



(24) 



1.4 X 1024 yrx/o,(Z,eo) 

Given the heavy neutrino mass Mn and the compositeness scale A, we only need 
to evaluate the nuclear matrix element Qfi to know the upper bound on the value 
of I/I imposed by neutrinoless double beta decay. 

The evaluation of the nuclear matrix elements was in the past regarded as 
the principal source of uncertainty in Oz//5/3 decay calculations, but the recent 



high-statistics measurement [^ of the allowed 2z//3/3 decay, a second order weak- 
interaction /5 decay, has shown that nuclear physics can provide a very good descrip- 
tion of these phenomena, giving high reliability to the constraints imposed by 0z//3/3 
decay on non-standard model parameters. 

Since we simply want to estimate the order of magnitude of the constraint in Eq. 
(24) we will evaluate the nuclear matrix element only approximatively. First of all 



the expression of the nuclear operator in Eq. (18) is simphfied making the following 
replacement 

' kV kl , ^ ' kV kl ^_^ -""A /ocN 

' kl ' kl ■^ 

The operator Q becomes then 

Q ^ .J^^rriARo) E r-{k)T^{l) (% - ^ak ■ ff] Fiv(x«) (26) 

where F^ = {1/x){Fb - 3Fa) = {l/A8)e-%x'^ - 3x - 3) with Fb and Fa given in 
Eq.(15). 

Since we are interested in deriving the lowest possible upper bound on |/| given 
by Eq. (24), let us find the maximum absolute value of the nuclear matrix element 
of the operator Q in Eq.(18): 

|%| < ^^imARo)\FN{x)\ \%\Mf\ + 1\Mgt\] (27) 

where MjT' = < f\J2k^i'r-{k)T_{l)\i > and M^r = < f\J2k^i'r-{k)r^{l)ak-ffi\i > are 
respectively the matrix elements of the Fermi and Gamow- Teller operators whose 
numerical values for the nuclear system under consideration are [^ |2^, Mp = 
and Mgt = —2.56. Inspection of the radial function Fn (for x > 0) shows that its 
maximum absolute value is attained at x = 0. In Eq. (27) we have evaluated F^ 
at a; = 2.28 {r^i = 0.5 fm). This value of r^^ corresponds to the typical internuclear 



distance at which short range nuclear correlations become important |T9|, so that 
the region x < 2.28 does not give contributions to the matrix element of the nuclear 
operator. We thus find 

IQfil < 0.6 X 10^ (28) 

which together with Eq. (24) gives the conservative upper bound on |/| shown in 
Fig. 2 as a function of M^ for A = 1 TeV. 

In particular, we see that the choice |/| ~ 1 is compatible with bounds imposed 
by experimental limits on neutrinoless double beta decay rates. We emphasize that 
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our bound on |/| is conservative, and an exact evaluation of the nuclear matrix 
element will give an even higher lower bound. 

We also note that Eq. (24) can alternatively be used to give a lower bound on A 
as a function of M^ (assuming |/| = 1 ). This is shown in Fig. 3 where we can see 
that the lower bound on the compositeness scale coming from Oz//5/5 decays is rather 
weak: A > 0.23 TeV at M^ = 1 TeV. In table I we summarize our bounds for some 
values of the excited Majorana neutrino mass. We remark that, as opposed to the 
case of bounds coming from the direct search of excited particles, our constraints 
on A and \f\ do not depend on any assumptions regarding the branching ratios for 
the decays of the heavy particle. 

To obtain more stringent bounds, we need to improve on the measurements of 
0^(3(3 half-life. However, our bounds c.f. Eq. (24) on (|/| or A ) depend on the 
experimental T1/2 lower limit only weakly (oc T^/2 ) so that to obtain an order of 
magnitude more stringent bound we need to push higher, by a factor of 10^, the 
lower bound on Ti/2- 

We should bear in mind, however, that the simple observation of a few OujSf] 
decay events, while unmistakably proving lepton number violation and the existence 
of Majorana neutrals, will not be enough to uncover the originating mechanism 
(including the one discussed here). In order to disentangle the various models, single 
electron spectra will be needed, which would require high statistics experiments and 
additional theoretical work. 
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TABLE CAPTIONS 

[Table I] Lower bounds on A with |/| = 1, and upper bounds on |/| with A = 1 TeV, 
for different values of the heavy neutrino mass Mn ■ 



FIGURE CAPTIONS 

[Fig. 1] Schematic illustration of neutrinoless double beta decay Oz//3/3 via the exchange 
of a Majorana neutrino. 

[Fig. 2] Conservative upper bound on |/| versus the heavy Majorana neutrino mass 
Mjv, at A = 1 TeV. 

[Fig. 3] Conservative lower bound on A versus the heavy Majorana neutrino mass Mjv 
for I/I = 1. 
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Figure 3 



